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ABSTRACT 

In  this  report  the  equations  of  motion  for  the  six 
degrees  of  freedom  of  an  aircraft  were  solved  by  an  analog 
computer  circuit.   These  solutions  were  incorporated  with 
the  instrumentation  of  a  fixed  base  flight  simulator  to 
provide  a  variable  stability  aircraft  demonstration  device 
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I.   INTRODUCTION 

This  report  is  the  fourth  in  a  series  of  projects  to 
convert  a  "Jet  Propelled  Aircraft  Instrument  Flying 
Trainer,  Type  C-ll"  manufactured  by  Link  Aviation,  Inc., 
Binghamton,  New  York,  to  a  variable  stability  fixed  base 
aircraft  simulator.   The  objective  being  the  provision  of 
a  dynamic  real  world  demonstration  device  for  use  by  the 
Aeronautics  Department  at  the  Naval  Postgraduate  School. 

The  first  report  [Ref.  1]  considered  the  feasibility 
of  converting  this  C-ll  for  use  as  a  variable  stability 
simulator.   It  was  concluded  that  the  computer  system 
available  could  be  made  compatible  with  the  simulator  if 

S/->Tnp      ovfoTiciiro     TTi<->/-1-i-P-i(~'3+--i(-inc      i.r<=>  -r<=»     m  n  A  p>      to      "f"V»  <=>      T.TTllc 

The  second  report  [Ref.  2]  concerned  the  physical 
alterations  to  the  trainer  and  the  replacement  circuits 
necessary  for  simulator- -computer  compatibility.   The 
overall  dimensions  of  the  C-ll  were  reduced  and  the  A-C 
computer  components  were  removed  and  replaced  by  D-C 
servo  systems. 

The  third  report  [Ref.  3]  integrated  the  longitudinal 
equations  of  motion,  as  solved  on  the  analog  computer,  with 
the  simulator  controls  and  instrument  displays.   The  con- 
clusions of  this  third  report  stated  that  the  quality  of 
the  flight  characteristics  of  the  trainer  were  not  satis- 
factory if  the  simulator  was  to  be  useful  as  a  demonstration 
device . 


The  goals  of  the  present  fourth  project  were  twofold. 
The  first  objective  was  to  improve  the  flying  qualities  of 
the  simulator  itself.   This  entailed  improving  or  modifying 
the  instrumentation  of  the  simulator  so  that  it  would  faith- 
fully represent  the  computed  solutions.   The  second  objec- 
tive was  to  integrate  the  full  six  degrees  of  motion  into 
the  computer  solution  and  hence  provide  a  complete  flight 
envelope . 


II.   THE  COMPUTER -SIMULATOR  ARRANGEMENT 

The  hybrid  computer  system  available  for  use  with 
the  simulator  consists  of  a  ComCor  CI-5000  Hybrid  Analog 
Computer,  an  SDS-9300  Digital  Computer,  two  Adage  AGT-10 
graphics  display  units  and  various  input  output  devices, 
(Figure  1) .   The  computer  systems  are  operated  and  main- 
tained by  the  Electrical  Engineering  Department  and  the 
cockpit  simulator  by  the  Aeronautics  Department  at  the 
Naval  Postgraduate  School. 

The  SDS-9300  (Scientific  Data  Systems  Corp,  El 
Segundo,  California)  has  a  32k  24-bit  core  memory,  two 
million  6-bit  character  drum,  two  magnetic  tape  units, 
line  printer  output,  and  multiple  input  peripherals. 
Its  application  to  the  real  time  solutions  required  for 
Link-computer  interaction  is  enhanced  by  three  features: 
a  priority  interrupt  structure  and  external-word  parallel 
input/output  channel  to  the  CPU,  and  a  second  port  to 
memory  which  allows  memory  access  at  a  500  Khz  word  rate 
without  degrading  central  processor  operation. 

The  ComCor  CI-5000  (ComCor  Inc.,  Anaheim,  California) 
as  presently  installed  has  24  Summing  Amplifiers,  28  Com- 
bination Amplifiers  (Summer  or  Integrator),  10  Multipliers, 
1  Resolver,  48  Servo  Set  Potentiometers,  32  Manual  Set  Po- 
tentiometers and  other  common  analog  components,  as  well 
as  a  limited  parallel  digital  logic  computer  with  "and" 
gates,  flip  flops,  counters,  and  multiple  clocks.   The 
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CI-5000  has  a  serial  character-oriented  adder  and  control 
system  specifically  designed  for  digital  computer  control. 

The  CI-5000  was  the  only  computer  used  directly  in  the 
problem  solution  for  this  project,  although  the  SDS-9300 
was  used  to  set  the  servo  potentiometers.   Meaningful  in- 
corporation of  the  graphics  display  unit  into  a  hybrid 
solution  is  recommended  as  a  future  project. 

Trunklines  are  hardwired  from  the  simulator  to  the 
CI-5000  and  may  be  patched  into  the  problem  solution  from 
terminals  on  the  analog  board.   Each  of  the  control  inputs 
is  taken  from  a  tapped  potentiometer  in  the  Link  and  sent 
to  the  computer  as  a  voltage.   And  in  turn  the  voltage 
outputs  from  the  computer  solutions  are  sent  to  the  Link 
■\kjQ tq    thev  a1"*^  tT*°ns^"0T*med  ^~  o  visual  re°dr,utc  ori  ^"hQ  slid ■" 
propriate  cockpit  instruments. 

The  system  originally  required  at  least  two  people, 
in  two  different  rooms,  to  operate  it.   One  was  needed  to 
start  and  stop  the  computer  and  the  other  to  fly  the  simu- 
lator.  This  inconvenience  was  overcome  by  using  a  logic 
circuit  (Figure  2)  that  allowed  the  computer  to  be  started 
and  stopped  from  a  switch  in  the  cockpit. 
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Figure  2.   Fly-Reset  Logic  Circuit 
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III.   AIRCRAFT  AND  AERODYNAMIC  EQUATIONS  OF  MOTION 

The  stability  and  body  axes  systems  used  to  describe 
the  force  and  moment  equations  were  determined  from  the 
data  and  information  available  prior  to  this  project. 
However,  a  brief  description  of  these  systems  will  be 
given  for  easy  reference. 

The  body  axes  system  is  located  at  the  center  of 
gravity  of  the  aircraft  and  fixed  with  respect  to  it. 
The  mutually  perpendicular  axes  have  the  conventional 
orientation  with  x  being  the  longitudinal  direction,  y 
the  lateral,  and  z  the  vertical.   The  determination  of 
whether  the  longitudinal  axis  parallels  the  fuselage  axis, 

Wiri  rr      r*  *"»  r\  7*H  r\  ir       t  nrnc  t       1    t  -n  o       to       n>*n'i'f">*iT*'\r       ort^l      t»to  q       r~  V*  *-^  <-  r>  -r>       +•% 

parallel  the  fuselage  reference  line. 

The  stability  axes  system  is  also  located  at  the 
center  of  gravity  of  the  aircraft  and  fixed  with  respect 
to  it.   However,  the  x  axis  is  chosen  so  that  it  is  posi- 
tive in  the  direction  of  motion  of  the  aircraft  with  re- 
spect to  a  steady,  symmetric,  reference  flight  condition. 

A  graphical  illustration  of  the  axis  system  is  shown 
in  Figure  3.  For  this  report,  the  x,  and  x  axes  will  be 
assumed  parallel  for  the  chosen  reference  flight  condition 

The  reference  angles  are  the  classical  Euler  Angles 
and  are  shown  in  Figure  3.   In  terms  of  the  body  axes  the 
angular  rates  are  given  in  Figure  4. 
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;.  yb 


zb   ze 


Figure  3.   Euler  Angles. 


<$>h    =   pb    +    tan    9b    (qb    sin    <f>b    +    rb    cos    <|>b) 
9b    =   qb    cos    ^    -    rb    sin   ^ 
%    sin    *b    +    rb    cos    *b 


* 


cos 


Figure  4.   Angular  Rates. 

The  equations  of  motion,  as  given  in  Reference  4,  are 
listed  in  Figures  5  and  6.   All  moment  equations  are  written 
with  respect  to  the  body  axes  and  all  force  equations  with 
respect  to  the  stability  axes. 

Lb  =  XxPb  +  <W  qb  rb  "  Jxz  (Pb%  +  V 


Mb  =  ly%    +    CIX"IZ)  P 


b  rb  "  Jxz  (rb  "  Pb} 


nk  =  Ilr-  +  (I  -I  )  p,  q,  -I   fp,  -  r,q,  ) 
b     z  b    *■  y   x^  FbHb.  xz^Fb    bHb^ 

Figure  5.   Body  Axes  Moment  Equations. 
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X  =  m  [u   -  v  (r,  cos  a  -  p,  sin  a) 
s      L  s    s v  b         *d      J 

+  g(sin  8,  cos  a  -  cos  8b  cos  <J),  sin  a)] 

-  T  cos  (a+aT) 

Yg   =   m    [vg    +   ug    (rb    cos    a    -   pb   sin   a)    -    g   cos    6b   sin   4>b] 

Zg    =    m    [vs(pb    cos    a   +    rb    sin   a)    -    ug (qb    -    a) 

-  g(cos    6b    cos    <j>b    cos    a   +    sin    6b    sin   a)] 

-  T   sin    (a+aT) . 

Figure  6.   Stability  Axes  Force  Equations. 

As  noted  under  the  computer  description,  only  one 
resolver  is  currently  available  on  the  CI-5000;  therefore, 
small  angle  approximations  arc  essentially  required  to 
handle  the  axis  rotation  problem.   Further  simplifications 
based  on  the  small  angle  limitations  are  listed  in  Figure 
7.   (Reference  4) 

1)  a,  3,  9  are  small  [cos  a  ~  1,  sin  a  ~  a] 

2)  Products  and  squares  among  a,  3,  p,  q,  r  are 
negligible 

3)  ug  :  V   tan  3  :  B    -^  -  3 

s 

4)  Thrust  Component   T  sin  (a+aT)  negligible. 

Figure  7.   Simplifications. 

Substitution  of  these  approximations  into  the  force 
and  moment  equations,  and  considering  the  Euler  angle  rates, 
leads  to  the  simplified  equations  listed  in  Figure  8. 
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MOMENTS : 


Lb    •    Jxz      - 

—   =  pb  "  —  rb 

X  X 


Mb 


-y  =  % 


b         xz  • 

—  =  rb  "  —  pb 

z         z 


FORCES: 


v  -  -|  -  gCO-a  cos  *b)  ♦  I 

v(3+rb-pba)  =  —  -  g  sin  ^ 

Z 
v("'qb)  =  "m  +  g  cos  ^b 


EULER  ANGLE  RATES: 

*b  =  Pb  +  °b  ^b 

0b  =  qb  cos  4b  -  rb  sin  ^ 

^b  =  qb  sin  *b  +  r  cos  *b 
HEIGHT: 

h  =  v(6b-a  cos  <J>,  -  3  sin  cf),  ) 
Av  =  /  v  dt 

Figure  8.   Simplified  Equations 


The  aerodynamic  equations  of  motion  for  the  T-33  air 
craft  are  given  in  Reference  4  and  listed  in  Figure  9. 
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AERODYNAMIC  MOMENT  EQUATIONS  -  BODY  AXES 

3L  3L 

L,/I      =    LQ    3    +   -~     oS+'Lp   +  L      r   +   -«£     ar 
bx          3               3a                   pr  r  3a 


+    LfiA   *A   +    \   6R 


VXz    =    V+V+i     ap    +    Nr   r   +    N        *A 


3N- 

6A 
+   Nr      6n    +   — s —  a   <5. 
6R     R  8a  A 


M,  /I      =    M      a    +   M"    a   +    M      q,     +    Mr      6^ 
b'    y  a  a  q   nb  6p      E 


AERODYNAMIC    FORCE    EQUATIONS    -    STABILITY  AXES 


X      =    -    a      STCM      +    CI      a   +    C'         a2] 
s  -o       -    D  u  D    2 

o  a  a 


q0  S   1-  [C    Av  +  C'  aAv  ♦  C'    a2Av] 
\   °  '    o        a        a 


Ys  =  mv  [YB  6  ♦  Y6r  «r] 


o      a        5- 

b 


-  qB  S   -   [C.   Av  +  CT   aAv] 
0   1  vo     Lo       La 


Figure  9.   Aerodynamic  Equations. 
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IV.   ANALOG  SOLUTIONS 

The  physical  data  for  the  constants  and  derivatives 
of  the  T-33  aircraft,  as  given  in  Reference  4,  were  inserted 
in  the  listed  equations.   The  scaled  analog  computer  cir- 
cuit solutions  were  then  determined  as  shown  in  Figures 
10  through  13.   The  potentiometer  values  are  collected 
and  listed  in  Table  I. 

An  interface  network  is  shown  in  Figure  14  with  the 
potentiometer  values  included  in  Table  I.   These  added 
circuits  were  necessary  to  make  the  scaled  computer 
variables  compatible  with  the  separately  scaled  simulator 
inputs  and  control  outputs . 

dinal  equations,  as  developed  in  Reference  4,  were  solved 
on  the  digital  and  analog  computers.   Both  solutions  pro- 
duced the  unacceptable  aircraft  response  shown  in  Figure 
15.   Initially,  the  same  results  were  obtained  for  this 
study.   To  get  the  response  shown  in  Figure  16,  the  coef- 
ficient of  drag  (Cn  )  in  X  was  multiplied  by  ten  to  damp 

JJg        s 

the  long  period  oscillation,  M  was  multiplied  by  one  hun- 

H 

dred  to  damp  the  short  period  oscillation,  M  was  multiplied 
by  ten  to  increase  the  frequency  of  the  short  period  oscil- 
lation, and  M.   was  multiplied  by  ten  to  maintain  a  response 
to  a  change  in  the  elevator  deflection  (6p) . 

Although  it  appears  that  gross  variations  have  been 
introduced  in  to  certain  stability  parameters  in  order  to 
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Figure    10.      Analog   Circuits. 
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Figure    11.      Analog   Circuits 
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Figure  12.   Analog  Circuits 
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Figure  13.   Analog  Circuits. 
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Mr      421 


Av      ^-v     10  rs51 

— 5 f  J "^  Indicated  Airspeed 


*         A25      10  r4i 


50 


O10    isy  x 
^  Indicated  Vertical    Speed 


:  A22      i.r*4 


Oiu  rvL 
P> Indicated  Angle   of  Attack 


424                  63 
2a  /^\ 10    [X Indicated   Gyro    Pitch 


400  52 


"7  A 


o— "-> 


t (     ) is — i   •>  maicatea   byro   £anK 


416        in   .36    .205 

2v— O— IH>-0- 


Indicated   Rate    of  Turn 


Ah  ,s53 

5Q-  r^  Indicated  Altimeter 


406       013 
rpra Q       [p> Indicated  Tachometer 


Figure  14.   Interface  Circuits. 
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SETTING 


POT 
NUMBER 


SETTING 


000 
001 
002 
003 
004 
005 
006 
007 
010 
Oil 
012 
013 
014 
015 
016 
017 
020 
021 
022 
023 
024 
025 
026 
027 


1329 
3220 
0644 
5192 
0519 
0100 
9660 
0100 

3220 
5000 
0519 
3220 
6100 
1433 
1137 
2381 
1219 
0298 
2161 
4704 
1084 
0750 
,1277 


030 
031 
032 
033 
034 
035 
036 
037 
040 
041 
042 
043 
044 
045 
046 
047 
050 
051 
052 
053 
054 
055 
056 
057 


.1776 
.1229 
.0462 
.0102 
.1604 
.1319 
.2820 
.0858 
.1170 
.0287 
.0135 
.0549 
.1207 
.0113 
.0313 
.0172 
.2701 

.0132 
.0454 
.9720 
.0397 
.0500 
.6375 


400 
401 
402 
403 
404 
405 
406 
407 
410 
411 
412 
413 
414 
415 
416 
417 
420 
421 
422 
423 
424 
425 
427 
437 


.1666 
.9999 


6A(.9000) 


.0466 
.2168 


6  (.0800) 


.0361 


5E(.7200) 


5500 

h 

h 

5600 
3333 

1111 
1650 

*A 


TABLE  I.   Potentiometer  Settings 
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Figure  15.   Initial  Response. 
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Figure  16.   Modified  Response 
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obtain  a  stable  solution,  discrete  digital  solutions  to 
the  equations  of  motion  indicate  that  the  increased  coef- 
ficients are  nearer  to  predicted  values  than  are  the  co- 
efficients as  listed  in  the  reference.   The  order  of 
magnitude  change  in  the  drag  coefficient,  for  example, 
causes  the  drag  value  to  fall  nearly  on  the  curve  of  drag 
coefficient  versus  angle  of  attack,  while  the  values  de- 
rived from  the  data  and  formulae  of  the  reference  do  not 
approach  the  curve. 

The  lateral  equations  also  proved  to  be  unusable 
initially.   Once  an  angle  of  bank  (<j>)  was  assumed,  the  yaw 
angle  (3)  would  increase  till  computer  overload.   The  solu- 
tion was  made  stable  by  multiplying  NR ,  by  ninety  to  get  an 
acceptable  frequency  for  the  Dutch  Roll  mode. 

The  changes  made  were  arbitrary,  but  necessary  to  get 
a  useable  solution  from  the  data  available.   The  circuits, 
as  modified,  were  considered  to  represent  the  aircraft  in 
its  normal  flight  mode.   Demonstrations  that  changed  the 
stability  characteristics  of  the  aircraft  used  this  mode  as 
the  reference  condition. 
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V.   FLIGHT  SIMULATOR  MODIFICATIONS 

The  C-ll  Instrument  flight  simulator  was  not  orig- 
inally designed  to  represent  any  particular  aircraft. 
The  primary  use  was  for  instrument  instruction,  not  as 
a  variable  stability  demonstrator.   Consequently,  no 
provisions  were  made  for  selectively  varying  chosen  per- 
formance parameters;  hence,  the  previous  projects  pro- 
ceeded to  reconfigure  the  C-ll  to  convert  it  to  a  variable 
stability  simulator. 

With  one  exception,  all  the  flight  instruments  had 
been  converted  to  a  dc-ac  servo  system.   In  these  systems, 
a  dc  motor  was  driven  by  a  differential  amplifier  that 
compared  the  input  signal  to  a  feedback  signal  from  a 
follow  up  potentiometer.   Motor  movement  was  sensed 
through  a  gear  train  by  an  AC  syncronous  servo  system 
that  positioned  the  cockpit  instrument.   A  block  diagram 
representation  is  shown  in  Figure  17. 
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Figure  17.   Block  Diagram  Instrument  Drive  System 


26 


There  were  two  basic  problems  with  the  system  as 
designed  and  installed;  it  did  not  respond  to  small  in- 
puts and  once  it  did  respond,  it  would  not  re- center  to 
the  same  position  each  time  when  the  signal  was  removed. 
The  two  problems  are  directly  related  in  that  the  dc 
drive  motors  need  a  finite  voltage  to  respond  to  a  sig- 
nal, be  it  input  or  feedback.   For  the  motors  available, 
this  voltage  was  0.5v — £>4v.   With  the  cost  of  better  dc 
motors  prohibitive,  it  was  decided  to  replace  as  many  in- 
struments as  possible  with  direct  current  meter  movements 
and  use  the  best  of  the  available  motors  for  the  remaining 
instruments . 

Consequently  the  airspeed,  angle  of  attack,  vertical 
speed  indicator,  and  acce] erometer  servo  systems  were  re- 
placed with  dc  meter  movements.   The  indicator  instruments 
were  removed  and  the  dials  that  the  pilot  sees  were  mounted 
on  the  new  movements  and  replaced  in  the  cockpit.   Thus, 
the  pilot  still  sees  the  conventional  aircraft  instrumen- 
tation.  To  provide  for  calibration,  the  output  of  a  trunk- 
line  from  the  computer  connected  to  a  potentiometer  which 
acts  as  a  voltage  divider  and  supplies  an  appropriate  vol- 
tage to  the  indicator.   (Figure  18.)   A  resistance  was 
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Figure  18.   D.C.  Meter  Indicator  System. 
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placed  in  series  with  the  potentiometer  to  prevent  inad- 
vertently burning  out  the  meter  during  the  calibration 
process . 

The  altimeter  circuit  was  modified  from  the  original 
design,  Figure  19,  by  replacing  the  single  turn  follow  up 
potentiometer  with  a  three  turn  potentiometer  and  then 
increasing  the  feedback  gain  to  the  amplifier.   This  ar- 
rangement (elimination  of  R: ,  R2 ,  and  C)  provided  a  good 
centering  characteristic;  and  also  improved  the  response 
of  the  system  as  shown  in  the  before  and  after  response 
of  the  follow  up  potentiometer  movement.   (Figure  20.) 
The  input  resistor  was  also  reduced  to  increase  the  allow- 
able altitude  variation  from  6000  to  10,000  feet. 

The  gyro  bank  system  was  treated  similiarly  to  the 
altimeter.   The  response  was  improved  by  changing  the 
motor  and  increasing  the  gain  of  the  feedback  from  the  fol- 
low up  potentiometer.   Although  the  response  is  not  com- 
pletely smooth,  it  is  satisfactory  for  demonstration  pur- 
poses . 

The  gyro  pitch  system  was  left  as  an  AC  servo  system 
due  to  the  size  limitations  on  the  motor  (it  mounts  in- 
side the  gymbaled  ball).   The  poor  response  of  the  orig- 
inal design  was  traced  to  the  omission  of  a  feedback  path 
for  the  operational  amplifier  and  improper  phase  relation- 
ships in  the  chopper-amplifier  circuit  (Figure  21) .   By 
inserting  the  feedback  path  and  optimizing  capacitor 
values,  the  problem  was  corrected  and  the  circuit  responds 
smoothly  and  accurately. 
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The  tachometer  and  gyro  bank  systems  were  treated 
similarly  to  the  altimeter.   The  responses  were  improved 
by  changing  motors  and  increasing  the  gain  of  the  feedback 
from  the  follow  up  potentiometers.   Though  the  responses 
are  not  completely  smooth,  they  are  satisfactory  for  mean- 
ingful demonstration  purposes. 


32 


VI.   CONCLUSIONS 

The  stated  goals  of  this  project  were  to  improve 
the  simulator  response  and  to  incorporate  the  six  degrees 
of  freedom  in  the  computer  solution.   Both  objectives  were 
accomplished  but  with  varying  degrees  of  success. 

The  instrumentation  of  the  simulator  now  responds 
accurately  and  smoothly  to  the  computer  output.   The  con- 
trol stick  has  a  small  dead  space  in  the  centered  position, 
though  this  does  not  appear  to  significantly  effect  the 
pilot's  response.   The  operation  of  the  simulator  is  con- 
sidered quite  acceptable  for  variable  stability  demonstra- 
tions . 

The  solution  of  the  equations  of  motion  with  the  data 
given  in  Reference  4  are  not  considered  acceptable.   The 
reason  for  the  arbitrary  changes  that  were  needed  to 
utilize  the  solutions  are  not  understood  and  make  com- 
parison between  the  simulator  and  the  aircraft  itself 
meaningless.   However,  changing  stability  parameters  did 
change  the  responses  of  the  simulator,  though  it  did  not 
appear  that  all  of  the  reactions  were  in  proportion  to  the 
changes  made. 

Consequently,  it  is  recommended  that  for  future  work 
with  the  simulator,  new  aircraft  data  be  obtained  and  that 
the  aerodynamic  equations  be  solved  in  accordance  with 
that  data.   If  the  original  solutions  would  correspond  to 
the  actual  aircraft  response,  more  confidence  could  be 
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placed  in  the  varied  responses  due  to  changes  in  the  air- 
craft stability  parameters. 

The  ease  with  which  aircraft  parameters  could  be 
changed  was  demonstrated.   This,  combined  with  the  ability 
to  program  a  sequence  of  changes  on  the  digital  computer 
warrants  further  effort  for  this  project. 
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